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in ■ Abstract 

Qh! We examine the measured rates for the decays D K^*^£u, B K^*^'tp^'^ 

and B K*^ in a number of scenarios, in the framework of the heavy quark 
pH I effective theory. We attempt to find a scenario in which ah of these decays 

are described by a single set of form factors. Once such a scenario is found, 
^ \ we make predictions for the rare decays B K^*k+e-. While we find that 

H ' many scenarios can provide adequate descriptions of all the data, somewhat 

surprisingly, we observe that two popular choices of form factors, namely 
monopolar forms and exponential forms, exhibit some shortcomings, espe- 
cially when confronted with polarization observables. We predict Br{B^ 
K0^+^-) = 6.4 ± 1.0 X 10"'^ and Br{B° i?*V+Ai") = 3.8 ± 1.3 x lO^^. 
We also make predictions for polarization observables in these decays. 
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I. INTRODUCTION 



The decays of heavy hadrons have recently received much attention in the hterature 
p|^l3|. From the experimental standpoint, these decays allow access to some of the funda- 



mental parameters of the standard model, such as the elements of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. Questions of CP violation, heavy-flavor oscillations and many 
others have added to this interest. 

From the theoretical standpoint, these processes, and the heavy hadrons themselves, 
allow various quark models of QCD, as well effective theories, to be tested. In particular, 
the heavy quark effective theory (HQET) has both been tested by experimental observations, 
and has played a major role in the extraction of \Vcb\ from experimental data ||T3[| . 

Much of the success of HQET has been in the treatment of decays from one heavy flavor 
to another, namely b —>■ c transitions. The effective theory is more limited in scope when 
applied to heavy-to-light transitions, such as c — ^> s or 6 — > m. Neverthless, as we will outline 
in a later section, the scaling behavior of the form factors that describe various weak decays 
can be deduced WM. This, in principle, allows the form factors for 6 — s transitions to be 
inferred from those for c — s. 

In this article, we assume the validity of the heavy quark symmetry and examine the 
decays D K^*^iu, B K'^*^ip^'\ and the recently measured B K*'j in a number 
of scenarios. In particular, we seek a scenario in which all of these decays are adequately 
described by a single set of form factors. A number of authors have performed similar 
analyses ||T5|-[T9|, using the decays D K^*Hv and B K^*^^'\ with varying results. 

Once we find a scenario that is satisfactory for the decays mentioned above, we examine, 
briefly, the decays B K^*H^i^. In this way, we hope to make reliable estimates of the 
absolute decay rates for these processes. In a later article, we will consider more details of 
these decays, such as forward-backward asymmetries. These decays, as well as the decay 
B K*'y, are particularly interesting as the short distance operators responsible arise first 
at the one-loop-level, and are therefore sensitive to new physics beyond the standard model 
p0|^3|. However, in order for any effects due to such new physics to be clearly identified, 
the long distance contributions that arise in the hadronic matrix elements must be well 
understood. As witness to this, we point out that the inclusive rate b —>■ sj is reasonably 
well understood, while the exclusive rate B — > K*^ has been predicted to be anywhere from 
2% to 40% of this inclusive rate Pl . 



In the case of the nonleptonic and the rare decays, there arises the crucial issue of the form 
factors describing the B K^*^ transition. Such form factors may be estimated in various 
models, from QCD sum rules, or by applying the scaling relations predicted by HQET to 
the form factors for the corresponding D — > K^*^ transition. The method suggested by 
HQET, while clearly model- independent, is itself somewhat problematic, as one must take 
into account the fact that predicting B —>■ K^*^'ip, for instance, requires that these relations 
be carefully extrapolated beyond the kinematic range accessible in the D K^*^ transitions. 
This is because the maximum in the D K^*^ transitions is 1.95 (0.95) GeV^, while the 
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appropriate to the B — > K^*^t/j transition is 9.6 GeV^. 

The HQET symmetry predictions relate form factors at the same value of the kinematic 
variable v-v' (or v-p'), not g^, where v and v' are the four- velocities of the parent and daughter 
hadron, respectively. In this case, the extrapolation is from v -v' = 2.1 m D ^ Kiu to v -v' 
= 3.6 in — > Kijj. The corresponding numbers for the K* decays are 1.3 and 2.0. The 
question of extrapolation also applies to the rare decays, and is particularly important for 
the decay B — > K*^, for which — 0, but v • v' — 3.0. 

For the nonleptonic decays, a second issue is that of the factorization approximation, 
which is commonly used to calculate the hadronic matrix elements required. This approx- 
imation is not very well founded in general theoretically, yet it appears to work well phe- 
nomenologically in the B decays where it has been tested. Nevertheless, application of the 
form factors of some model or effective theory to the decay B Kip, in conjunction with 
the factorization approximation, serves to probe both issues, and may fail due cither to 
inadequate choice of form factors, failure of the factorization approximation, or both. 

The rest of this article is organized as follows. In the next section we describe briefly the 
effective Hamiltonians and the appropriate hadronic matrix elements for the processes of 
interest. In section III we use HQET to obtain relations among the form factors of interest. 
In section IV we discuss our fltting procedure, and present the results for the decays that 
we are considering. In section V we discuss possible limitations of our results, and suggest 
questions that may be of interest to both theorists and experimentalists. 



II. DECAY PROCESSES 



A. Semileptonic Decays 

Of the three processes we discuss, the semileptonic decays are perhaps the simplest to 
treat theoretically. The effective Hamiltonian for these decays is 

Tieff = ^VcsSl^ (1 - (1 - 75) J^i. (1) 

The hadronic matrix elements for the decays D —>■ K^*^^ are 

{K{p') \s^,c\ D{p)) = U{p + p')^ + - p'),, 

{K{p')\s^^^^c\D{p)) =0, 

{K*{p\ e) \s^,c\ D{p)) = ige,,^pe*^{p + p'Tiv - p'f. 

{K*{p',€)\s-f^-f5c\D{p)) = f€*^ + a+e* ■p{p + p')^ + a.€* -pip-p')^. (2) 

These decays are thus described in terms of six independent, a priori unknown form factors. 
The terms in /_ and a_ are unimportant when the lepton mass is ignored, since 

(p - P')/7''(1 - 75)^^£ = {k, + h)^h^{l - 75)^/^ = mdl'^il - l^W (3) 
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In experimental analyses, the form factors /+, a+, / and g are usually assumed to have 
the form 

M^') = (4) 

where = {p — p'Y ■, and is a mass, usually taken to be that of the nearest resonance 
with the appropriate quantum numbers. To date, only the mass appropriate to /+ has been 
measured, and its value measured by CLEO is 2.0 ± 0.22 GeV [^. The masses appropriate 
to a+, / and g are assumed to be 2.5 GeV, 2.5 GeV and 2.1 GeV respectively. The /j(0) 
have the values 0.66 ± 0.03, -0.14 ± 0.03, 0.40 ± 0.07 and 1.55 ± 0.11, for /+, a+, g and /, 
respectively p6 . 



B. Nonleptonic Decays 



Neglecting penguin contributions, the effective Hamiltonian for the nonleptonic decays 
of interest here is 



Kff = [Ci{m,) (C7^(l - 75)&) {si'{l - 75)c) 

+ C2{m,) (S7^(l - 75)&) (CT^I - 75)c)] , 



(5) 



where 



Ci(mb) 

C2(mb) 



-6/23 



a 



12/23' 



as{niw) ) ' \as{niw) ) 

-6/23 / / N \ 12/23' 



a 



a 



as{mw) J 



(6) 



Here, mw is the mass of the W boson, and Q;s(/i) is the running coupling of QCD. This 
effective Hamiltonian mediates two classes of B nonleptonic decays. The first class contains 
a D in the final state: B DX where X may be Dg, D*. The second class contains a light 
meson in the final state: B 'i^'' Y where Y is now a charmonium state. 

To evaluate the matrix elements of the effective Hamiltonian we employ the factorization 
assumption. By Fierz rearrangement we rewrite the effective Hamiltonian in a form which 
is suitable for use with this assumption. Both terms of the effective Hamiltonian contribute, 
in general, but for the decays in which we are interested, only the second term is of interest, 
and it may be written 



cS 



C2{mb) + —Climb) 



X (s7^(l - 75)6) (c7^(l - 75)c) 



(7) 



where Nc is the number of colors. 
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At this point, it has become customary to replace the couphng coefficient by a phe- 
nomenological constant 02, whose absolute value is measured to be approximately 0.24. The 
sign of 02 is not important for our discussion at this point. It will become important if long 
distance contributions to the dileptonic rare decays are included. 

For the decay B — ^ K^*\j/-ip, we therefore write, after using factorization, 

J/^K'-*^ \n,s\ b) = ^V^bKla^ {J/^P |c7m(1 - 75)c| 0) 



X 



K(-*^\sYil-l5)b\B). (8) 



The hadronic matrix elements (^i^^*^ |s7'^(l — 75)6! are analogs of those of the previous 
subsection, so that the form factors required for these matrix elements may, in principle, be 
obtained from the 'semileptonic' decays of B mesons into kaons. Such decays do not take 
place in the standard model, but we can invoke heavy quark symmetries to relate the form 
factors needed to those for the semileptonic decays of D mesons to kaons. 
The remaining matrix element is 

{J/^p{p^, e^) IcYC^ 75)c| 0) = Um^e'^. (9) 

The decay constant can be obtained from the leptonic width of the appropriate charmo- 
nium vector resonance as 



U = f^-^-^'^^ (10) 

where we have ignored lepton masses. In this way, we find /^=0.382 GeV, and /^/=0.302 
GeV. 



C. Rare Decays 



In the standard model, the effective Hamiltonian for the decay 6 ^ 57 is pD|,|22 

Gf e 



K 



eff 



V*VtkCrimt)sa^, K (1 + Ts) + (1 - 75)] 



(11) 

where F'^'^ is the electromagnetic field strength tensor, and the term in may be safely 
ignored. For the decay b —>■ si^i' , the corresponding effective Hamiltonian is 

Gp a 

72 477 



.rUb 



2^^C7(m5)sc^^,g^(l + 75)6^7^^ 



+ C9(m,)s7^ (1 - 75) bi^i + Cio(mfe)s7^ (1 - 75) blj^j.i 



(12) 



The Wilson coefficients Ci are as in the article by Buras et al. [^. For the discussion at 
hand, we have ignored contributions that arise from closed qq loops, although these may be 
easily included. 
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The only new hadronic matrix elements that arise in the rare decays of interest are 



{K{p') \7sa^,b\ Bip)) 



IS 



[p + p')^ {p - p')^ -{p + p')^ {p - p')^ 
{K\p\ e) \m,,b\ B{p)) = e,.^p [g^e*'^ {p + p'f + g.e*'^ {p - p'f 

+he* ■ pip + p')°' {p — p')^ 

introducing four new form factors. Due to the relation 



(13) 



(14) 



we can easily relate the matrix elements above to those in which the current is 5(7^,^756. 
Experimentally, nothing is known about the form factors s, g± nor h. 



III. HQET AND FORM FACTORS 

Using the Dirac matrix representation of heavy mesons, one may treat heavy-to-light 
transitions using the same trace formalism that has been applied to heavy-to-heavy transi- 



tions p8| , |29| . In the effective theory, a D meson traveling with velocity v is represented as 

2ig 



DM 



-75 = Md{v), 



(15) 



with an identical representation for a B meson. The meson states of the effective theory are 
normalized so that 



{Viv')\Viv)) = 2vo6'' (p-p'). 
The states of QCD and HQET are therefore related by 



(16) 



(17) 



In all that follows, we will represent the states of full QCD as \D{v)), and the states of 
HQET as \V{v)). 

Let us consider transitions between such a heavy meson [D meson) and a light meson 
(kaon) of spin J, through a generic flavor changing current. In the effective theory, the 
matrix element of interest is 

K^-^\p,v) \sThi''^\v{v)) = Tt{ETMd{v)}. (18) 

r is an arbitrary combination of Dirac 7 matrices, and 77 is the fully symmetric, traceless, 
transverse, J-index tensor that represents the polarization of the state with spin J. The 
matrix S must be the most general that can be constructed from the kinematic variables 
available, and Dirac 7 matrices. The most general form for this is PD 
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,75 



(19) 



The ^j's are uncalculable, nonperturbative functions of the kinematic variable v ■ p, and the 
1 (75) is present if the resonance K^'^^ has natural (unnatural) parity. 

From this point on, let us limit the discussion to only two of the kaon resonances, namely 
the ground-state pseudoscalar kaon itself, and its vector counterpart, K*. The above then 
leads to 



Kip) 
K*ip,e) 



V{v] 
Viv 



Tr {[(^3 + ^^^4) t*-v+^* (^5 + m] TMd{v)} . 



(20) 



At this point, let us emphasize that the form factors are independent of the form of F, 
and so are valid for decays through the left handed current (F = 7^(1 — 75)), as well as for 
rare decays (F = cr^^(l ±75)). These form factors are also independent of the mass of the 
heavy quark, and are therefore universal functions. Thus, they are valid for D K decays, 
as well as ioi B K decays. This independence of the quark mass allows us to deduce, 
in a relatively straightforward manner, the scaling behavior of the usual form factors that 
describe these transitions \T^\. We illustrate this by examining one matrix element in detail. 

Consider 



{Kip) |s7^c| D{po)) = f+{po+p)f, + /-(Po -p)a 



= ^/mB[K{p) 
From these equations, one finds that 



V{v)) = ^Tt I {^1+^^2)151, 



1 + 



(21) 



6 = ^(/. + /- 



6 



(/- - /+) 



(22) 



Since the do not scale with the mass of the heavy quark (or meson), it is trivial to deduce 
that 



(23) 



For the other transitions of interest, the form factors are as defined in the previous 
section, and the relationships between these and the are 



J, = (/,+/.) 
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if- - /+) = -VmBs, 



^3 
^4 
^5 



m 



3/2 
D 



(a+ + a_) 



{2g — a+ + a_) = m'Jy'^h 
(/ + 2mDt; ■ pg) -- 



(9, 



9-] 



= ^/mD9 



1 



{9- - 9+) 



These expressions yield the additional scaling behavior 



(24) 



-3/2 

D 5 



J. 1/2 -1/2 



-1/2 



m 



-1/2 



, h ^ m^, 
1/2 



D 



3/2 
9+- 9- 



m 



D 



1/2 



Eqns. (|T^) and (^OD , and consequently eqns. (|2Tl - p^ , contain all of the leading order mn 



dependence in the form factors, and are valid irrespective of the mass of the strange quark. 
If the strange quark could be treated as heavy, then we could think of the as arising from 
an infinite sum of terms in the l/m^ expansion. The leading order forms (in l/m^) are also 
contained in these expressions. These expressions are also valid in the limit of a very light 
strange quark. 

Isgur and Wise [|1^] have used the scaling of /+ + /_ (~ ) to say that this combi- 

nation of form factors vanishes (at leading order in 1/m/)), and suggest that one can write 
/_ = — /_|_ + (9(l/m£)). Implicit in this argument is the assumption that the strange quark is 
very light. In our formalism, this amounts to setting ^1 to zero, and we would automatically 
lose the full scope of our predictions. This is because we could then never recover the limit 
of a heavy strange quark, for in this limit, ^1 = — i/m^^/2, where ^ is the usual Isgur- Wise 
function. 

The strange quark is such that it may be treated as either heavy or light. We believe that 
neither the full heavy s limit (^1 = — i/rn^^/2), nor the limit of a very light s quark (^1 0) 
is completely satisfactory. We assume neither limit in our analysis, and therefore make full 
use of the predictions of the heavy quark effective theory, which are valid independent of the 
mass of the strange quark. This means that we retain the form factor ^1 in our discussion 
and treat it as a completely independent form factor, tying it to neither of the two limits 
discussed. 

At the risk of boring the overly patient reader to tears, and perhaps even to death, we list 
one more set of relationships among form factors, this time writing the usual form factors 
in terms of the ^j. The relations are 
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1 

1 



(^1 - 

(6 + mD^2) 



+ = ^ [6 + (6 - ^4) 



m 



D 

1 



3/2 

D 



9 



6. 



(25) 



For the corresponding transitions with a b quark (and a B meson) in the initial state, all 
factors of mo above must be replaced by m^. Using this, rearrangement of eqns. (^) yields 



f ■ p) 

V ■ p) 

V ■ p) 

V ■ p) 

V ■ p) 

V ■ p) 







n 





















f-iv-p) 



1 + 



TUB J 

tub) 



+ f-{v-p) 
+ f+{v-p) 



rriB 
rriB 



1 



- 1 



ffnB\ 
KnipJ 
frriD 



1/2 



r{v-p), 



1/2 



1 /m^y/^ 

2 v^y 

1 

2 KniB 



9 (v-p), 



) 



n 



{v ■p)[l + 
(vp) (1 + 



rriB 
mp 

rriB 



+ a^{v-p) 



+ a'livp) 



mp 

rriB 
mp 



- 1 



(26) 



mB 

where is the form factor appropriate to the D ^ K transition, while is the form 
factor appropriate to the B K transition, and quantities on the left-hand-sides of eqns. 
(p6|) are evaluated at the same values of f ■ j9 as those on the right-hand-sides. Omitted from 
each of eqn. (^) is a QCD scaling factor, discussed below. 

Eqns. ( p6D illustrate two effects, namely the scaling of form factors in going from the D 
system to the B system, as well as the mixing of a+ with a_, and /+ with /_. The effect of 
this mixing is very important in going from D transitions to B transitions, as it introduces 
form factors that have not yet been measured experimentally, or to which experiments are 
not yet sensitive, namely /_ and a_. In the rates for the semileptonic decays D — > K^*^^, 
terms dependent on /_ and a_ are proportional to the mass of the lepton, and thus play a 



miniscule role, except near 



Such terms may also be significant if the polarization 



of the charged lepton is measured. 

We close this section with a brief discussion of radiative corrections to the currents and 
matrix elements that we have discussed. In the limit of a heavy b quark, the full current of 
QCD is replaced by ^ 
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sTb 



(27) 



This arises from integrating out the b quark, and matching the resulting effective theory 
onto full QCD at the scale m^, at one loop level. At the scale m^, we must also integrate 
out the c quark, but there is also the effect due to running between nib and nic. The net 
effect of this is that the form factors C,i appropriate to the b ^ s transitions are related to 
those for the c — > s transitions by 



_6_ 

' 25 



(28) 



The forms of the matrix elements that we discuss above are valid in the limit of infinitely 
heavy b and c quarks. For quarks of finite mass, there are clearly going to be corrections to 
the relations we have obtained. In other words, new form factors that appear first at order 
l/nic and 1/mf, will begin to make contributions. It is expected that such contributions 
will become more significant away from the 'non-recoil' point, or for f ■ p > m^{*). This is 
particularly important for the B — > K^*^ decays, as v-p can become very large. Nevertheless, 
we will apply the relations we have found through all of the available phase space. It is our 
hope that by doing this, we will at least be able to indicate the suitability of HQET for 
these transitions. However, since we fit the form factors rather than attempt to calculate 
them in some model, some of these higher order effects may, in fact, have been included. 



IV. RESULTS AND DISCUSSION 
A. Data 

All of the results we describe are for decays with (or K^*) in the final state. The 
treatment of charged kaons would be identical, and we believe that our results in these chan- 
nels would be similar in quality to those we obtain for neutral kaons. Before describing the 
results of our fits, we must comment on how we treat the available data, particularly in the 
case of the semileptonic decays. Very few of the experimental collaborations have extracted 
acceptance-corrected distributions [E3,R2]. Instead, the form factors are all assumed to be 



of the monopole form p6[, and the parameters are then extracted from the Monte Carlo 
simulations, with acceptances and efficiencies folded in. 

Because of this, we proceed in the inverse sense to generate some 'simulated data'. We 
use the published monopole parameters for the form factors to generate dV /dq^ spectra for 
the semileptonic decays, using the published uncertainties in the monopole parameters to 
generate uncertainties in the simulated data. In general, these errors are correlated, but we 
ignore this correlation. 

The simulated data generated in this way are completely smooth. We introduce an 
'anti-smoothing' by smearing the simulated data with a pseudo-randomly generated gaussian 
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distribution of mean zero and standard deviation determined by the errors in the unsmeared 
simulated data. It is this smeared simulated data that we use for fitting. For the decays 
K^*i^i', we also include the ratios F^/r and F+/F_ in the fit, where F± are as 
defined in PDG [^. In addition, we must point out that the measured decay rates for 

— ^ K^^*H~^u are somewhat smaller than those for — > K~^*^i^u, while the published 
form factor parametrizations are for the average over charge states. To account for this, we 
rescale the values of the /i(0) to correspond to the smaller rates for neutral kaons. It is these 
rescaled values that are cited in section |II A| , and that we use to generate the simulated data. 

For the nonleptonic decays, we fit to the PDG averaged decay widths for K^ijj 
and K'^*ip. In the case of the latter, we also include the ratio F^/F. For this ratio, 

we take the averaged value of 0.78 ± 0.073 as calculated by Gourdin et al. |T5|,|T6[. Masses 
and lifetimes of mesons are all taken from PDG |2^, and we use Vtb = 0.9988, Vts = 0.03, 

= 0.9738, Kfe = 0.041, mfe = 4.9 GeV, m,=1.5 GeV, mt=m GeV. 

It is worth mentioning that the experimental choice of monopole form factors may be 
inappropriate, particularly for /. In the limit of a heavy strange quark, one finds that 
/ oc {1 + V ■ v') S,{v ■ v') , where ^ is the Isgur-Wise function. If ^ is assumed to be monopolar 
in g^, then simple pole dependence for / is inappropriate. Even in the case of a light strange 
quark, we find that / oc ^5 + w ■ p^g? again indicating that the dependence on the kinematic 
variable is not simply a monopole form. However, for the range of accessible in D ^ K* 
decays, the decay rate is not sensitive to the kinematic dependence. The effect on the decay 
rate for the nonleptonic and rare decays being considered here, however, are quite significant. 

In our fitting, we have separated the decays containing a K meson in the final state from 
those containing a K* meson in the final state. Thus, for instance, we do not include data 
on ratios of rates like r{D ^ K*ih')/r{D Kiu). Our reason is that such ratios introduce 
correlations between the parameters of the form factors for the two sets of decays. 



B. Heavy-s Limit 

One approximation used recently in the literature has been to treat the strange quark 
as heavy [p3| , p^ , so that the decays of interest can be treated in the heavy-to-heavy limit. 



In this limit, the form factors of eqn. ( pOD may be written 

^1 — ^5 — ^ ^; 

^2 = ^6 = „ , C, 

6 = ^4 = 0, (29) 

where ^ is the Isgur-Wise function for heavy-to heavy transitions, and in this limit, rriK = 
iriK*- In particular, this limit means that ^5 = ttik^g- 

We have used this form in fitting the data, and have obtained reasonable fits to the 
differential decay rates in the semileptonic decays, as well as to the total decay rates in 
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the nonleptonic decays. Polarization ratios, however, are poorly reproduced. In the case of 
the nonleptonic decay B K*ip, the ratio Tl/T depends only on kinematics, and has a 
value of 0.43, independent of the form chosen for ^. The experimental value is 0.78 ± 0.07. 
In addition, the ratio r_|_/r_ in D ^ K*ii> always has a value of about 0.4, significantly 
different from the experimental value of 0.16±0.04, and largely independent of the form 
chosen for ^. This indicates that the value of 0.43 obtained for the ratio Tl/T in B K*ip 
is not necessarily due to the breakdown of factorization in the heavy s limit, as this limit does 
not even provide an adequate description of all measurements in the semileptonic decays. 

Relaxing the strict heavy-s limit, but constraining the form factors to be near this limit, 
does not help much, as the polarization observables are still not well reproduced. The 
conclusion that we draw from this is that the heavy-s limit may give an acceptable description 
of unpolarized data, but may be dangerous when applied to polarization observables. 



C. General Features Of Results 



All of the results we present are obtained by performing four kinds of fits, namely (1) 
include the semileptonic decays D — > K^*^£u only; (1) include the semileptonic decays as well 
as the nonleptonic decays B K^*^ip; (3) include the semileptonic decays, the nonleptonic 
decays B K^*^^, and the nonleptonic decays B K^*^tp'] (4) include the measured decay 
rate for B K*'~f in the fit, as well as the three other decays. Clearly, in the case of the 
decays to K mesons, we need perform only three kinds of fits. In cases where a measured 
quantity is not included in the fit, we calculate that quantity using the fit parameters. 

We have explored two sets of parametrizations of the form factors. In the first scenario, 
which we refer to as the exponential scenario, each has one of the forms 



= ai exp [~hi {v -p- ni^i*))] = ai exp 

= ai exp -bi {v-p- mj^(,)f = ai exp 
^i = a^exp -bi {v-pf 



2m o 

bi. 



^max 



4m 



2 _ 2 
2 V ^max q 
D 



(30) 

(31) 
(32) 



while in the second scenario, which we call the multipolar scenario, the forms chosen are 

6 = ai{l + hv ■ pp , (33) 



with Hi = —2, —1, 0, 1. In the exponential scenario, eqn. (^) most closely corresponds to 
the forms that arise in some quark models, most notably that of Isgur and collaborators . 
However, in such models, the exponential of eqn. ( |30[ ) is usually multiplied by a polynomial 
in V ■ p, or equivalently, in g^. In the multipolar scenario, = —2 corresponds to a dipole 
form factor, while = — 1 represents a monopole. In any one fit, we do not choose all 



the form factors to have the same form. This means that, for instance, in the case of the 
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decays to K mesons, for which there are two form factors, the second scenario corresponds 
to sixteen different possible combinations of forms for and ^2- 

In each case, the Oj and bi are the free parameters to be varied in the fit. There are 
therefore twelve free parameters in each fit, to be compared with five extracted and three 
assumed parameters in the measured semileptonic decays. Since we have more free param- 
eters than are extracted from the experimental data, one might expect that it should be 
very easy to account for all of the data. In fact, the number of free parameters poses some 
problems, as it means that the problem is not very well constrained. One consequence of 
this is that there appear to be several local minima for any particular choice of form fac- 
tors. Nevertheless, there are some combinations of form factors that simply do not provide 
adequate descriptions of the data, despite the large number of free parameters. 

When only the semileptonic decays are included in the fit, we find that almost any com- 
bination of forms for the form factors leads to reasonable results. The few combinations that 
do not provide good descriptions fail only in their description of the polarization observables. 

When we use the exponential forms, we find that we are not able to obtain an adequate 
description of all of the data simultaneously. In particular, when we include the nonleptonic 
decays in the fit, the prediction for the rare decay B K*'-f is significantly different from 
the measured rate. In some cases, however, we find that if we fit the semileptonic decays 
alone, omitting all of the nonleptonic decays, the prediction for the rare decays is of the right 
order of magnitude. One possible conclusion here is that factorization is not applicable to 
these nonleptonic decays, or that there are significant non-factorizable contributions to the 
amplitude. 

In contrast with the exponentials of the first scenario, many combinations of the 'multi- 
polar' forms of the second scenario lead to good descriptions of all the data simultaneously. 
One outstanding feature of all of our results in this scenario (which also exists in the expo- 
nential scenario, but to a lesser extent) can be easily understood by examining eqn. (^5|), 
where it is seen that is present in all of the form factors that describe the semileptonic 
decays D —>■ K* . It is therefore not at all surprising that our results are most sensitive to 
this form factor. Invariably, we have found that the best fits occur for linear mv-p {i.e., 
ne = 1), independent of the forms chosen for the other form factors. Furthermore, the slope 
parameter feg is almost always negative, with values lying between -0.4 and -0.65 GeV~^. 
The only positive values for fog occur when only the semileptonic decays are included. 

The fact that is so well constrained in our fits means that / and g are also quite 
well constrained. Since these are the only two form factors that are needed for the decay 
B K*^, it is no surprise that we find little variation in the predictions for this decay rate 
as we vary the parametrizations of ^3, .^4 and .^5, provided that the choice of parametrization 
of ^6 is unchanged. 

For other forms of we find that at most two of the semileptonic, nonleptonic or rare 
decays are well accomodated. For instance, if we choose a monopole form for then in 
addition to the semileptonic decays, we find that we can accomodate either the nonleptonic 
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decays, the rare decay, but not both. In addition, if we choose all form factors to be 
monopolar, we fail to find adequate descriptions of the polarization observables, particularly 
for the ratio Ti/T measured in the nonleptonic decay. 

Our results for the form factors are comparable with those of other authors. Gourdin et 
al. |[TB| have found that a number of scenarios, including monopolar form factors, are unable 



to describe the nonleptonic measurements. Aleksan et al. JT^ have found that softening the 
scahng relations allows an adequate description, while in a second article, Gourdin et al. fl^ 



have found that allowing the form factor / to decrease linearly with allows an adequate 
description of the data. We find that / is quadratic in g^, but the absolute value at = 
for the B K* transitions is larger than the absolute value at = ^max) keeping with 
the results of [|16[, and quite different from pole models. 



D. Parameters And Form Factors 

For each scenario, we have selected a set of fits that we consider to be representative. 
The parameters for these fits are displayed in table |I| for the exponential forms, and in table 
|l| for the multipolar forms. In the first scenario, the results we have selected correspond to 
^2 and ^3 as in eqn. (^), and ^4 as in eqn. (pi]), and ^5 and as in eqn. ( |3^). In the 
second scenario, the Ui have the values Ui = (0, —1, 1, —2, 1, 1). 
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TABLES 

TABLE L Values of the parameters that result from four different fits, for the exponential 
scenario. In this table, Fit 1 means that only D K^*H V is included in the fit; Fit 2 means 
D K'^*ki^ and B K^*\j/^ are included; Fit 3 means D K^*kv, B K^*'^J/'4) and 
B K^*^i;' are included; Fit 4 means D KHu, B K*J/ip, B K*ip' and B K*-y are all 



included, and applies only to decays with K*^s in the final state. 



Parameter 


Fit 1 


Fit 2 


Fit 3 


Fit 4 


ai 


2.497 


0.671 


0.667 




12 


1.503 


-0.265 


-0.266 




13 


10.0 


8.825 


1.745 


1.742 


a4 


4.521 


5.137 


1.259 


1.257 


as 


9.996 


9.998 


9.990 


9.726 


aa 


1.044 


1.037 


1.080 


1.074 


bi 


3.353 


1.979 


1.984 




b2 


4.710 


8.0xl0~''^ 


3.1x10"^ 




bs 


3.511 


0.238 


7.6x10^5 


3.0 xlO"'* 


bi 


10.0 


1.447 


1.076 


1.076 


bs 


5.914 


6.730 


6.520 


6.458 


be 


0.585 


0.563 


0.585 


0.581 



TABLE II. Values of the parameters that result from four different fits, for the multipolar 



scenario. The columns are as in table 



Parameter 


Fit 1 


Fit 2 


Fit 3 


Fit 4 


ai 


3.358 


0.144 


0.158 




(12 


2.929 


5.875 


5.951 




13 


-10.0 


0.162 


0.144 


0.145 


a4 


-10.0 


3.154 


2.569 


2.365 


15 


1.399 


0.288 


0.346 


0.451 


16 


0.055 


1.094 


1.085 


0.922 


bi 










b2 


0.290 


10.0 


10.0 




bi 


-0.507 


10.0 


10.0 


10.0 


bi 


1.134 


0.468 


0.388 


0.354 


bs 


-1.012 


-1.042 


-1.030 


-0.907 


be 


10.0 


-0.439 


-0.435 


-0.392 
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To give some sort of meaning to these parameters, we display in tables |T| and |I^, the 
values of the form factors at = 0, as well as their logarithmic derivatives at the same 
point. In these tables, the theoretical error that we quote, as well as those that we quote for 
the rest of our results, are estimates only, and are obtained by using the covariance matrix 
that results from the fit. 

For the monopole form chosen by experimentalists, the logarithmic derivative is 



1 



m. 



2 ' 



(34) 



where mj is the polar mass. In this way, we can compare our form factors at = 0, and the 
corresponding slope parameters, with the experimental values. It is gratifying to find that 
for /, g and /+, the values we have obtained are, for the most part, quite consistent with 
the experimental values, for both the exponential and multipolar scenarios. In addition, 
the logarithmic derivative of a+ shows some variation as we go from fit to fit, and there are 
sizable variations in the values for /_, a_ and their logarithmic derivatives. These variations 
are not surprising, as there are no experimental constraints on these quantities. Although 
the values we obtain for these quantities may lie outside the accepted domain suggested by 
models, it is nevertheless quite satisfying to note that the values do not change much as we 
go from Fit 2 to Fit 3 to Fit 4, particularly for the multipolar scenario. It is also interesting 
that the prefered slope parameter for a+ is negative. 

TABLE III. Values of the form factors and their logarithmic derivatives at = 0, for the 



Quantity 


Experiment 


Fit 1 


Fit 2 


Fit 3 


Fit 4 


/+{0) 


0.66±0.03 


0.59±0.04 


0.66±0.03 


0.66±0.03 




4(0) 

/+(0) 


0.23 ± 0.09 


0.80±0.19 


0.24±0.07 


0.24±0.07 




/-(O) 




0.97±0.48 


-0.06±0.12 


-0.07±0.12 




fL(0) 

f-W 




0.98±0.17 


-2.50±4.61 


- -2.32±3.96 




/(O) 


1.55±0.11 


1.53±0.10 


1.55±0.07 


1.57±0.11 


1.57±0.08 


f'(0) 
f(0) 


0.16 


0.15±0.12 


0.12±0.05 


0.14±0.03 


0.14±0.05 


9(0) 


0.40±0.07 


0.35±0.03 


0.36±0.02 


0.37±0.03 


0.37±0.02 


9'(0) 
9(0) 


0.23 


0.36±0.12 


0.35±0.05 


0.36±0.02 


0.36±0.05 


a+(0) 


-0.14±0.03 


-0.23±0.05 


-0.20±0.05 


-0.19±0.05 


-0.19±0.04 


a'^(O) 
a+(0) 


0.16 


-3.17±0.66 


-1.72±0.76 


0.03±0.11 


0.03±0.44 


a_(0) 




-2.99±0.36 


-6.33±1.98 


-1.18±0.05 


- -1.18±0.28 


a'_(0) 
a_{0) 




1.25±0.15 


0.12±0.08 


- -0.05±0.01 


-0.05±0.11 
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TABLE IV. Values of the form factors and their logarithmic derivatives at = 0, for the 
multipolar scenario. The columns are as in table |. 



Quantity 


Experiment 


Fit 1 


Fit 2 


Fit 3 


Fit 4 




U.ODitU.Uo 


U.D4diU. Uu 


n fiQ-un no 


U.DzitU.U4 




/+(0) 


0.23 ±0.09 


0.29zb0.17 


0.29ib0.07 


0.29±0.08 


- 






-O.ODlti.oo 


n Q /I _l_n Q /I 
-U.o4itU.o4 


n Qfi-Ln Q/1 




/-{(>) 




0.03±0.03 


0.21±0.04 


0.21±0.04 




/(O) 


1.55±0.11 


1.55±0.14 


1.55±0.08 


1.53±0.18 


1.54±0.07 


f'(0) 
f(0) 


0.16 


0.04±0.12 


0.14±0.06 


0.17±0.20 


0.15±0.05 


9(0) 


0.40±0.07 


0.51±0.05 


0.40±0.03 


0.40±0.05 


0.37±0.03 


9'(0) 
9(0) 


0.23 


0.22±0.01 


0.24±0.03 


0.23±0.05 


0.19±0.02 


a+(0) 


-0.14±0.03 


-0.24±0.05 


-0.21±0.04 


-0.20±0.08 


-0.20±0.03 


a'^(O) 
a+(0) 


0.16 


-1.14±0.33 


-1.11±0.30 


-0.93±0.71 


-0.98±0.27 


a_(0) 




3.54±0.33 


-1.37±0.25 


-1.21±0.61 


- -1.21±0.22 


a'_{0) 
a_{0) 




0.22±0.02 


-0.08±0.01 


- -0.09±0.11 


-0.09±0.01 



The form factors for the D K^*^ transitions that result from these fits are shown in 
fig. |l| for the exponential forms, and in fig. |^ for the multipolar forms. In each of these 
figures, we indicate with arrows the maximum possible in the corresponding D K^*^ 
decay. The corresponding form factors for the B K^*'^ transitions are shown in figs. |^ 
and § respectively. In these latter figures, the vertical dashed lines indicate the range of 
that correspond to the D K^*^ decays: the values oi v ■ p within the dashed lines of figs 
and ^ are the same as those that lie to the left of the arrows in figs. |l| and ^ 
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FIGURES 




FIG. 1. Form factors for the D K^*^ transitions that result from fits, using exponential 
forms. The thick solid curve is the experimentally extracted form factor, while the thick dashed 
curves show the range that results from the experimental uncertainties in the parametrizations of 
the form factors. In each graph, Fit 1 means that only D K^*^(,v is included in the fit; Fit 2 
means D K^*ku and B K^*U/ijj are included; Fit 3 means D K^*ku, B K^*U/tJj 
and B K'^*)ij;' are included; Fit 4 means D K*£u, B K*J/ip, B K*ij;' and B K*^ 
are all included, and applies only to a+, / and g. On the scale shown, the curves for Fit 4 are 
indistinguishable from those of Fit 3. The arrows indicate the maximum accessible in the 
semileptonic decays D K^*ku. 
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FIG. 2. Form factors for the D K^*^ transitions that result from fits, using 'multipolar' 
forms. The key is as in fig. |^. 

The effect of the mixing mentioned near the end of the previous section is seen in the 
curves for and a^. In the case of /+ all of the curves for the D form factors are very 
close to each other, and are all quite similar to the monopole form over the entire range of 
physically accessible g^. Application of eqn. (^61) , or more truthfully, of eqn. (pO|), leads to 
curves for seen in figs. || and |. We point out that the scaling effect due to the coefficient 
of the term is only about 15%, so that the very different forms for seen in the figures 
must be attributed to the mixing with f^. 

In examining the curves of this figure, one must remember to compare the form factors 
near the kinematic end points. This means that near = g^ax^ — {''^d — — 1-95 

GeV^ should be compared with near = {mB—TriK)'^ = 24.8 GeV^. It is very interesting 
to note that the form factors for the B K^*^ transitions that result from fitting to 
the semileptonic decays alone are usually very different from those that result when the 
nonleptonic decays are included in the fit. 

By examining the graphs of figs. |l] and |], as well as the numbers of tables |Tl| and |^ it is 
clear that the form factors that we have obtained are quite consistent with the experimentally 
extracted ones, in the range of physically accessible q^. Outside of this range, however, all of 
the form factors we obtain are markedly different from the experimental forms. This clearly 
has very important consequences for the nonleptonic and rare decays. 
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FIG. 3. Form factors for the B K^*^ transitions that result from fits, using exponential forms. 
In each graph, Fit 1 means that only D K^*^ii' is included in the fit; Fit 2 means 
and B K^*'^J/tlj are included; Fit 3 means D K^*ku, B K^*^J/ip and B K'^*^' are 
included; Fit 4 means D K*iv, B ^ K*J/tp, B ^ K*tp' and B K*^ are all included, and 
applies only to a+, / and g. The region between the vertical dashed lines is the range of for 
which information is available from the D K^*kv semileptonic decays. 
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FIG. 4. Form factors for the B K^*^ transitions that result from fits, using multipolar forms. 
The key is as in fig. ^ 

E. Results And Discussion 

As can be seen from the numbers in tables |V] and |V|, all fits to the semileptonic decays 
are reasonable. The quality of the fit with respect to the nonleptonic decays, the results 
of which are displayed in tables VI I| and [VI I H , is quite different, however. In the case of 



Fit 1, the nonleptonic decays are generally poorly described, while in the case of Fits 2, 
3 and 4, the theory does a reasonable job of describing all of the data. The differences 
between Fits 1 and 2 are shown most graphically in the form factors of figs. ^ and ^. These 
differences have very significant effects on the decay rates for the rare decays. In general, 
the differences in the form factors among Fits 2, 3 and 4 are much less striking. We also 
point out that the striking differences in form factors are achieved, for the most part, with 
very small adjustments to the intercepts and slope parameters, as seen in tables |TT| and [rV| . 

TABLE V. Results of fits for semileptonic decays, exponential forms. The columns are as in 
table | . 
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Quantity 


Experiment 


Fit 1 


Fit 2 


Fit 3 


Fit 4 


ro^Kt. (lO-i-* GcV) 


4.16 ±0.50 


4.26 ± 0.28 


3.91 ±0.24 


3.91 ±0.24 




T'd^K'I. (10-1'' GeV) 


2.98 ±0.25 


2.88 ± 0.12 


2.87 ±0.14 


2.76 ±0.12 


2.77 ±0.12 




1.23 ±0.13 


1.19 ± 0.12 


1.23 ±0.13 


1.06 ±0.08 


1.06 ±0.10 




0.16 ±0.04 


0.20 ±0.02 


0.18 ±0.02 


0.18 ±0.03 


0.19 ±0.02 



TABLE VI. Results of fits for semileptonic decays, multipolar forms. The columns are as in 



table | . 



Quantity 


Experiment 


Fit 1 


Fit 2 


Fit 3 


Fit 4 


Tn^Kt. (10-1-* GeV) 


4.16 ±0.50 


3.96 ± 0.28 


3.94 ±0.23 


3.95 ±0.25 




rD_K*te (10-1^ GcV) 


2.98 ±0.25 


2.88 ± 0.12 


2.87 ±0.14 


2.85 ±0.13 


2.88 ±0.10 




1.23 ±0.13 


1.21 ± 0.10 


1.22 ±0.09 


1.20 ±0.13 


1.23 ±0.09 


^{D^ K*tu) 


0.16 ±0.04 


0.16 ± 0.03 


0.16 ±0.03 


0.17 ±0.04 


0.20 ±0.03 



TABLE VIL Results of fits for nonleptonic decays, exponential forms. The columns are as in 



table | . 



Quantity 


Experiment 


Fit 1 


Fit 2 


Fit 3 


Fit 4 


T-B^Ki, (10-ifi GeV) 


3.29 ±0.95 


1.27 ±15.0 X 10-* 


3.29 ±0.95 


3.33 ±0.95 




TB-,K*i^ (10-16 GeV) 


6.93 ± 1.33 


1.61 ± 1.13 


7.06 ± 1.48 


6.75 ± 1.11 


6.77 ±1.30 




0.78 ± 0.073 


0.13 ± 0.21 


0.77 ±0.08 


0.78 ±0.05 


0.78 ±0.07 


Tb^k^- (lO-i'' GeV) 


< 3.5 


1.57 ± 10.79 X 10-=* 


1.07 ±0.25 


1.08 ±0.25 




^B-,K*i^' (10-16 GeV) 


6.14 ±3.95 


22.90 ±9.10 


5.75 ±2.64 


9.46 ± 0.53 


9.46 ±1.18 






0.07 ± 0.04 


0.96 ±0.02 


0.76 ±0.03 


0.76 ±0.05 



TABLE VIIL Results of fits for nonleptonic decays, multipolar forms. The columns are as in 



table B. 



Quantity 


Experiment 


Fit 1 


Fit 2 


Fit 3 


Fit 4 


Vb^K4, (10-16 GeV) 


3.29 ± 0.95 


77.02 ± 59.28 


3.79 ±0.51 


3.82 ± 0.65 




^B^K'i, (10-16 GeV) 


6.93 ± 1.33 


280.0 ± 52.0 


6.96 ± 1.30 


6.86 ± 1.34 


7.69 ±1.11 




0.78 ±0.073 


0.88 ± 0.03 


0.78 ±0.07 


0.78 ± 0.07 


0.71 ±0.05 


TB^Ki,- (10-16 GeV) 


< 3.5 


26.76 ± 17.10 


1.73 ±0.33 


1.75 ± 0.38 




TB^K'i,' (10-16 GeV) 


6.14 ±3.95 


26.00 ±4.50 


9.11 ±0.97 


8.33 ± 2.81 


8.49 ± 0.80 






0.69 ±0.07 


0.71 ±0.05 


0.70 ±0.08 


0.69 ±0.04 
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One comment on the ratio = F^^/r in i? ^ K*ip is worth making. For many fits, 
we obtain values for this ratio that are essentially unity. Gourdin et al. |T5| , p!6| state that, 
assuming factorization of the transition amplitude, this ratio has a maximum value of 0.833. 
They go on to point out that observation of a value of pi greater than this value would be 
indication of significant non-factorizable contribution to the transition amplitude. In their 
notation. 



PL 



a - hxf + 2(1 + c2y2) 
with a, h and c being determined by kinematics, and 



(35) 



X 



y 



(36) 



The form factors A2 and V are those defined by Bauer et al. I^B]. The error in their 
claim arises from ignoring the possibility that x may be large, so that pi — 1- In our 
notation, x ~ (m^ — mK*Ya+l f ■ Thus, we would claim that measurement of pi greater 
than 0.833 is not an indication that factorization is breaking down. Indeed, we have assumed 
factorization, but in fit 1 of table |V11I| , we find a value of pi = 0.88 > 0.83. 

In tables |X|and ^ we display our predictions for the rare decays B —>■ K*'~f, B — > Ki~^i~ 
and B K*i^i~ , obtained using the form factors from the four different fits. The lepton 
spectra are shown in figs. ^ and |^. As expected, the form factors from Fit 1, especially 
in the exponential scenario, lead to rates that are ruled out by experimental observations, 
particularly in the case of the decay B —>■ K*'~f. In the multipolar scenario, the predictions 
from Fits 2 and 3 are somewhat larger than the CLEO measurement of B ^ while 
the result from Fit 4 is consistent with the measured value. 

TABLE IX. Predictions for decay rates of rare processes, exponential scenario. The columns 
are as in table ffl 



Quantity 


Experiment 


Fit 1 


Fit 2 


Fit 3 


Fit 4 


Tb^K'j (10-" GeV) 


1.76 ±0.83 


0.01 ±0.03 


0.02 ±0.02 


0.01 ±0.01 


0.01 ±0.01 


^B^K^+^- (10-i« GeV) 


< 158.0 


0.30 ±0.72 


0.08 ±0.02 


0.08 ±0.02 








0.12 ±0.05 


0.13 ±0.03 


0.13 ±0.02 


0.13 ±0.02 






0.13 ±0.02 


16.3 ±12.1 


1.15 ±0.20 


1.15 ±1.60 


^B-,K'u+u- (10-'' GeV) 


< 10.1 


0.25 ±0.06 


16.5 ± 12.1 


1.27 ±0.19 


1.27 ± 1.60 
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TABLE X. Predictions for decay rates of rare processes, multipolar scenario. The columns are 



as in table 



Quantity 


Experiment 


Fit 1 


Fit 2 


Fit 3 


Fit 4 


re-.K'-, (10-" GeV) 


1.76 ±0.83 


77.90 ± 17.80 


6.17±3.31 


6.34 ±4.35 


2.21 ±0.72 


^B^Ku+^' (10-i« GeV) 


< 158.0 


6.33 ± 5.65 


0.27 ±0.05 


0.28 ±0.05 








5.05 ± 1.04 


0.33 ±0.13 


0.34 ±0.15 


0.22 ±0.03 






72.7 ± 14.5 


1.82 ±1.49 


1.42 ±1.25 


1.45 ±0.54 


^B-,K'^+^- (10-18 GeV) 


< 10.1 


77.7 ± 14.4 


2.15 ±1.49 


1.76 ±1.26 


1.67 ±0.55 
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Given the failure of the exponential scenario to explain the B K*'y data, one might 
be tempted to discard its predictions for the dileptonic decays. However, we see that the 
integrated rates are, for the most part, quite similar to those predicted in the multipolar 
scenario. The spectra that result from the two scenarios are very different, however, and the 
predictions for the relative amount of longitudinally polarized and transversely polarized K* 
produced are also somewhat different. The exponential forms predict Tt/Tl ^ 0.1, while 
in the multipolar scenario, the ratio ranges between 0.15 and 0.24. The moral here may 
be that the exponential forms may be adequate for predicting total rates, but not decay 
spectra, nor polarization observables. We limit our discussion to the multipolar scenario in 
what follows. 

10"" r . 1 . 1 n 10'" F ' 1 ' 1 ' 1 ' q 



10"" r 




q' (GeVl q- (GeV') 

FIG. 5. Differential decay rates for the processes B K^i^^^ and B K* ^'^ , in the 
exponential scenario. The graphs are, starting at the top left and moving clockwise: B Kfi~^fi~; 
B — > K*^'^fi~; B K* jj,'^ ^~ for longitudinally polarized i^*'s; B — > K*^'^^~ for transversely 
polarized K*^s. In each graph, the key is as in fig. ^. 
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FIG. 6. Differential decay rates for the processes B — > Kfi^iJ,~ and B — > K* ji^ ii~ , in the 
multipolar scenario. The graphs are, starting at the top left and moving clockwise: B — > Kn~^iJ,~; 
B K*iJ,^fj,^; B K* for longitudinally polarized K*^s] B — > K* jj,^ for transversely 
polarized K*^s. In each graph, the key is as in fig. ^. 

The predictions for the process B Ki~^i~ are two to three orders of magnitude smaller 
than present experimental upper limits, while those for B — > K*i~^i~ are smaller than the 
experimental limits by factors of four to six. We point out, however, that the calculated 
rates for these last two processes do not include possible contributions from charmonium 
resonances, which will certainly alter the shape of the lepton spectrum, and should also 
increase the total decay rate. An investigation of this effect will be left for a future article. 
However, in at least one experimental analysis, kinematic cuts are imposed on the total 
mass of the lepton pair, so that events that may arise from either of the first two vector 



charmonium resonances are excluded p4|. In any case, our results suggest that the exclusive 
dileptonic decay to the K* should be observed in the near future. 

The results that we have obtained here again illustrate that the fit to the present 
D — > K^*^ semileptonic spectra alone is inadequate for providing information on B ^ K* 
processes. Even when we include the nonleptonic decays, the predictions for different decay 
modes (particularly B K*i^£~ with longitudinally polarized K*'s) are sensitive to the 
nonleptonic modes we include in the fit. If any of these predictions are to be taken seriously, 
we would suggest that most attention be paid to the predictions of Fit 4 in the multipolar 
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scenario, as this is the only scenario that adequately describes all of the data available. 

One of the features of the predicted spectra are the minima in the differential decay rates. 
These minima are the result of zeroes in the respective helicity amplitudes, and the question 
of whether or not these zeroes do indeed exist, and of their exact locations, will have to 
await a i?-factory. However, long-distance effects, such as those that arise from charmonium 
resonances, or even from the charmonium continuum, will at least alter the positions of the 
zeroes, and may wash out the effect altogether. 



V. CONCLUSION 

We have used the scaling predictions of HQET, together with the most recent data 
on D ^ K^*Hu semileptonic decays, to extract the form factors that describe the D — *■ 
K^*^ and B K^*^ processes. The latter we have applied to other processes, namely the 
nonleptonic decays B K^*^ip and B K^*'>tp', as well as the rare decays B — >■ K*'y and 
B K^*H+i-. In the case of the nonleptonic decays, we have assumed factorization of the 
transition amplitude is valid. We have also performed simultaneous fits of the semileptonic, 
nonleptonic and rare processes, and have found that HQET, together with factorization, 
provide an adequate framework for describing the observations. Our predictions for the 
modes B — > K^*^i'^i~ suggest that these should be measurable in the next generation of 
experiments, and certainly at the proposed B factory. 

Perhaps the greatest shortcoming of our fit procedure lies in how we handle the data, 
or simulated data, for the semileptonic decays. Ideally, we should have attempted to fit 
our choices of form factors to the experimentally measured differential decay rates. As a 
second choice, our choices of form factors should have been input into the experimental 
Monte-Carlo programs to obtain the fit parameters. In any case, it is clear that the present 
data, particularly in the D — > K*lv mode, are inadequate to sufficiently constrain the form 
factors. In addition, the differences in form factors between Fit 1 and Fits 2, 3 and 4 are 
quite striking. 

The scenario that best describes all of the experimental data is the multipolar one, and in 
this scenario, we find that the universal form factor is linear in v-p. Using this scenario, we 
predict Br{B^ K^^i+^i-) = (6.4±1.0) x 10"^ and 5r(S° ^ K*°fi+fi-) = (3.8±1.3) x lO'^. 
These numbers are consistent with other model calculations We also predict Tt/^l in 
50 ^ K*^fj,+^- to be 0.15 ± 0.07. 

To fully constrain the predictions of HQET, information on the form factors a_ and /_ 
is needed from the semileptonic decays. Such information can only be obtained from high- 
precision measurements of the decay spectra at low values of g^, particularly for semileptonic 
decays to muons, as well as by measuring the polarization of the charged lepton, again 
preferably the muon. In addition, the precision and statistics in the spectrum must be 
improved so that the form factor parameters can be extracted from the data, particularly 
for the decays D — > K*iv. Perhaps the ideal experiment would be the equivalent of present 
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CLEO experiments, in which the machine is tuned to be a source of BB pairs, produced 
from the strong decays of the T(45'). For D decays, the equivalent would be to produce a 
copious number of ■0(3770) 's, which can be realized at the proposed tau-charm factory. 
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